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INTRODUCTION 



> 

Semiclassical approaches to finite Fermi systems provide a very efficient way of extracting 
the average behaviour of relevant physical quantities which characterize such systems. The 
most well known example is the celebrated Droplet Model and its extensions developed by 
Myers and Swiatecki which describe nicely the average behaviour of the nuclear masses 1]. 
Semiclassical techniques have also been applied to study the average behaviour of other 
properties such as inertias 2], charge radii sj, one- and two-body matrix elements^], etc. 
Our aim here is to present a new Thomas-Fermi (TF) theory, i.e. the % — > limit, for 
describing the average trends of the effect of pairing correlations in finite Fermi systems. 
Semiclassical approaches to the pairing problem can be of interest in scenarios like, e.g., 
cold atomic gases where the huge number of particles makes the full quantal calculation 
numerically very complicated. Also these calculations are useful if one is only interested in 
the average behaviour of the pairing gap, as is the case of the pairing term in the nuclear mass 
formula. The Local Density Approximation (LDA) is the standard semiclassical technique 
for dealing with the average behaviour of the pairing which was developed by Schuck and 
collaborators more than twenty years ago [10]. In LDA one considers the BCS equations in 
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infinite homogeneuous matter and replaces the Fermi momentum hp by its local version in 
terms of the density. The validity of LDA applied to the pairing problem is restricted on one 
hand to situations where the local Fermi wavelength 27r//cp(R) is small as compared with 
the distance where the mean field potential varies appreciably. In the case of a harmonic 
oscillator potential V(R) = mu 2 R 2 /2 this distance is the so-called oscillator length defined 



as I = ^Jh/muo. On the other hand, a second length scale introduced by pairing is the 
coherence length £ which measures the extension of the Cooper pairs. The validity of LDA 
in the pairing case also implies that the coherence length be smaller than the oscillator legth, 
i.e. £/Z < 1, which is usually equivalent to the condition A/hu) > 1 , where A is the gap 
in the single-particle spectrum. The condition £/Z < 1 is always violated in the outer tail 
of the surface because the LDA coherence length behaves as £ ~ A^ 1 and the gap vanishes 
in this region. In spite of these deficiences, integrated quantities as pairing energies may be 
quite accurate when considered on average [lol ]. 

In this contribution we present a novel TF theory for pairing which improves the LDA. 
This theory can be applied in the weak pairing regime where the chemical potential /i and 
the Fermi energy Ep have similar values and A//x << 1. This TF theory works, for the 
average, in the region A < hu, where LDA generally fails. 

An important point that will be discussed in this contribution is the possible quenching 
of the pairing gap when approaching the drip line. This fact is well documented in the 
recent nuclear physics literature [5 9 ] and we will see, using some examples, that this effect 
is a general feature when superfluid (superconducting) fermions come in a finite confining 
potential to an overflow situation. 

The contribution is organized as follows. The basic theory is presented in the second 
section. The main results are discussed in the third section. Our conclusions are laid out in 
the last section. 



II. BASIC THEORY 



It is well known that in the Hartree-Fock-Bogoliubov theory single-particle density matrix 
and the pairing tensor or anomalous density matrix are simultaneously diagonalized by the 
so-called canonical or natural basis, \n c ) As far as in this work we are only interested 
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in the weak coupling limit where the gap is small as compared with the Fermi energy, i.e. 
A//i < 1, one can replace with only small error the canonical basis by the basis of the normal, 
non-superfTuid mean field (HF) Hamiltonian. In this situation the gap equation reduces to 
its BCS approximation and can be written as 

A n = - ^ Vnn> T^hi W 
n' n 

where V nn > = (nn\v\n'n') is the matrix element of the interaction with \fi) the time reversed 
state of |n) and E n = [(e n — p) 2 + A 2 ] 1 / 2 the quasi-particle energies, with e n the diagonal 



11] written in the basis of the standard mean 



elements of the normal mean field Hamiltonian 
field hamiltonian, that is H\n) = e n \n). 

At equilibrium and for time reversal invariant systems canonical conjugation and time 
reversal operation are related by (r\n) = (n\r) =>- (ri^lrm) = (ri|p n |r2), where p n = \n){n\ 
is the density matrix corresponding to the state \n). Therefore the pairing matrix element 
can be writeen as: 

Km' = (nn\v\n'n') = J dr 1 dv 2 dr' 1 dr' 2 (r\\p n \r 1 )(r 1 r 2 \v\r' 1 r' 2 )(r' 2 \pn\^2)- (2) 

The density matrix p n fulfills the Schodinger equation 

(H - e n )p n = 0, (3) 

therefore we can write A n = Tr[Ap n ] and e n = Tr[Hp n ] and consequently the state depen- 
dence of the gap equation (CQ) is fully expressed through the density matrix p n . 

Performing the Wigner transform (WT) of Eq. (j3]) and taking into account that the WT 
of the product of two single-particle operators A and B equals, to lowest order in H, the 
c-number product of the corresponding WT's, i.e. A(R, p)B(R,, p), one easily obtains the 
h -»■ limit of Eq.© [ll | 

(if c ,-e)/ e (R,p) = 0, (4) 

2 

where H c i = 2m » ^ + V"(R) is the classical Hamiltonian which contains a local mean field 
potential V(R) and a position dependent effective mass m*(R) and / e (R, p) is the Wigner 
transform of p n . Equation (j3J) has to be read in the sense of distributions. Taking into 
account that x5(x) = one obtains the normalized distribution function 

f E (R, p) = -y^s5(E - H cl ) + 0(h 2 ), (5) 
9 (A) 



which corresponds to the Thomas-Fermi (TF) approximation of the normalized on-shell or 
spectral density matrix [4]. Its norm is equal to the level density (without spin-isospin 

degeneracy) : 

g TF (E) = — 1— f dRdp6(E - H cL ). (6) 



The semiclassical pairing matrix element can then be written as |4j: 



V ^ E>) = Iw^J ^/,(R,p)^(R',p')-(R,p;R',p'), (7) 

where v(R, p; R', p') is the double WT of < r 1 r2|f|r / 1 r / 2 >. For a local translationally 
invariant force this matrix element reduces to v (R, p; R', p') = 5(R — R')t>(p — p') with 
v (p — p') the Fourier transform of the force v(r — r') in coordinate space. 

The gap equation in the TF approximation is obtained by replacing in (JTJ p n and V nn > 
by their corresponding semiclassical counterparts Eqs. fl5]) and (J7|) respectively. In this way 
the TF gap equation reads 

ME') 



poo 

A(E)= / dE'g TF (E')V(E,E') 
Jo 



2 v /(E'-fiy + A*(E'y 

Eqs. ([H])-© can readily be solved for a given mean field and the chemical potential is fixed 
by the usual particle number condition. 

III. RESULTS 

As a realistic application of our TF theory we analyze the semiclassical pairing gaps 
as a function of mass number along the tin isotopic chain from 100 Sn to 132 Sn. To this 
end we use the D1S Gogny force [12| for both, mean field and pairing fields. The main 
ingredients for solving the semiclassical pairing equation (jHJ) are the on-shell density matrix 
/e(R, p) ©, which depends on the classical Hamiltonian H c[ that is determined by the 
effective mass m*(R) and the mean field V^R). These two quantities, namely m*(R) and 
V(R), are obtained through the Extended Thomas- Fermi (ETF) theory for finite-range non- 
relativistic interactions Q, Q]. Using these quantities as input, one obtains the level density 
([6]) and the pairing matrix element (j7]) which allow to solve the gap equation (jSJ) in our TF 
approximation. As explained in Refs. fl^ . Jjj], the ETF energy density functional can be 
transformed, inspired by the Kohn-Sham scheme, into a quantal functional from where the 
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FIG. 1: Left: Average pairing gap along the Sn isotopic chain. See text for details. Right: Average 
TF gaps at the Fermi energy as a function of the chemical potential in a WS potential computed 
in a box of radius R=25, 50, and 100 fm. 

quantal average gaps are obtained. It should be noted that within this approximation the 
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quantal functional associated to a finite- range effective interaction becomes local 
The quantal pairing gaps averaged with u 2 v 2 are depicted by circles in the left panel of Fig. 
1 and show the typical arch structure. In the same panel we also display the semiclassical 
TF gap at the Fermi energy by a thick solid line. We see that in this case the quantal arch 
structure completely disappears, as expected due the absence of shell effects in this case, 
and that the TF gaps decrease smoothly when the neutron number increases. As it has 
been discussed in Ref. [16], quantal effects, i.e. the arch structure over the shell, can be 
almost recovered by introducing some additional quantal fluctuations in the level density 
and retaining the TF pairing matrix elements ([7]) in the gap equation ([8]). The average gaps 
obtained in these conditions are displayed by diamonds in the left panel of Figure 2. We see 
that the arch structure for the tin isotopic chain is recovered and that the quantal gaps are 
predicted quite well in this way. 

An important feature of the average TF gaps is that they show a downward trend with 
increasing neutron number along a given isotopic chain as it can be appreciated in the left 
panel of Fig. 1. To investigate the behaviour of the TF gap at the Fermi energy approaching 
the drip line, we show in the right panel of Fig. 1 the TF gap as a function of the chemical 
potential employing a Woods-Saxon potential as mean field 17| . From this panel it can be 



seen that in the TF limit the gap vanishes just at the drip line when the chemical potential 
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FIG. 2: Left: Energy per baryon as a function of the average density in different WS cells corre- 
sponding to the inner crust of neutron stars. Right: Radial dependence of TF gap in the analyzed 
WS cells. End points indicate radius of WS cells. 



equals the depth of the single-particle potential. Actually we put the Woods Saxon potential 
in a large external container with radius R what can simulate a Wigner-Seitz (WS) cell as 
relevant in the inner crust of neutron stars, see below. We see that the gap becomes more and 
more suppressed at the edge of the Woods-Saxon potential for increasing values of the cell 
radius. Continuing filling the cell with a neutron gas, the gaps raise again. This quenching 
of the gap when approaching the drip line has also been studied by Hamamoto [^j analyzing 
the effective gap in weakly bound neutron levels in spherical and deformed nuclei, finding 
that the presence of a S1/2 component in the wavef unction of these levels strongly reduces 
the effective gap. Other gap components with I ^ show a less decreasing tendency near 
the drip line probably because of the centrifugal barrier which keeps the wave functions 
localised. This could mask the behaviour of the state dependent gaps in real nuclei near the 
neutron drip line. 

As mentioned already, another scenario where the quenching of the neutron gap appears 
is near the neutron drip line in the inner crust of neutron stars. This region is a crystal of 
nuclei embeded in a gas of free neutrons and electrons. The inner crust was described by 



Negele and Vautherin 



18] at HF level by means of the energy density functional method 



together with a spherical WS approach to deal with the crystal structure in an approximated 
way. The WS cell is electrically neutral and the ground state of the system of neutrons, 
protons and electrons is reached when they are in /3-equilibrium. We have performed a 
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similar calculation but at TF level and using the BCP energy density functional [19j. This 
functional consists of a bulk part provided by a microscopic calculation complemented by 
a phenomenological surface term. This functional with only four adjustable parameters 
reproduces nuclear binding energies and charge radii of finite nuclei with the same quality 
as obtained with the most performant effective forces. The ground state energy per baryon 
obtained in this way for average densities in the WS cells ranging from 2.79 x 10~ 4 to 
7.89 x 10~ 2 fm~ 3 are displayed in the left panel of Fig. 2 (black dots) in comparison with 
the Negele-Vautherin results (red diamonds) finding an excellent agreement between both 
calculations. Although pairing correlations in the inner crust of neutron stars are mainly 
driven by the free neutron gas, they have, however, a noticeable influence on the composition 
and pairing properties of the nuclear cluster inside the WS cell Q. On top of this TF 
calculation in the inner crust, we have also performed a TF pairing calculation in the studied 
WS cells. To this end, we have used the Gogny D1S force renormalized by a factor 0.85 to 
take into account that in the BCP functional the effective mass equals the physical one. A 
specially relevant quantity in this context is the radial dependence of the pairing gap A(R) 
which is obtained from the WT of the average quantal gap A av = Tr[Ak}/Tr[k] 20] which 
reads: 

A(R) = ^k)/(W A(R ' p)K(R ' p) (9) 

where k(R, p) = / dEg TF (E) K (E)f E (R, p) and A(R, p) = - / ^j»(p - p>(R, p'), 

We see in the right panel of Fig. 2 that A(R) takes a constant value in the outer part of 
the WS cell which corresponds to the gap of the free neutron gas 9J. The predicted behaviour 
for A(R) are in qualitative agreement with previous calculations ^6 9 ] . The nuclear cluster 
inside the WS cell disappears when the homogeneous phase is reached at an average density 
of about 0.08 fm -3 ( 982 Ge32). In this case the gap in the cell is practically the gap obtained in 
pure neutron matter at the same density. The density of the free neutron gas diminishes in 
approaching to drip configurations. In this situation the gap is strongly reduced not only in 
the gas but also inside of the nuclear cluster ( 200 Nb 4 i) and it even may disappear completely 
in the nucleus when the drip line is reached (see right panel of Fig. 1). Therefore, locally 
the TF A(R)'s are qualitatively different from what LDA would predict. 
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Conclusions and Outlook 

We have presented a TF theory for pairing in finite Fermi systems for weak coupling 
situations where A/ef << 1- This TF theory differs from the usual LDA. This essentially 
stems from the fact that we approximate the gap equation in configuration space and, thus, 
keep the size dependence of the matrix elements of the pairing force. This is not the case in 
LDA where the matrix elements of the force are always evaluated in plane wave basis. This 
semiclassical approach to pairing is only based on the usual validity criterion of Thomas- 
Fermi theory, namely that the Fermi wave length is smaller than the oscillator length. At no 
point the LDA condition that the coherence length must be smaller than the oscillator length 
enters the theory. Thus, the present TF approach yields for all pairing quantities the same 
quality as TF theory does for quantities in the normal fluid state. An interesting feature 
of our study is that the average gap breaks down going to the drip line. This unexpected 
result is confirmed by quantal calculations, though strongly masked by shell fluctuations. 
For systems with large numbers of particles the fluctuations should die out and, thus, the 
semiclassical behaviour prevail. Indeed preliminary results in a slab configuration show good 
agreement between quantal and TF gaps around the drip region. We also investigated in slab 
geometry the inverse scenario where the external potential gets in the upper part suddenly 
strongly constricted rather than widened. Very preliminary results show that the gaps now 
become much enhanced where before they were suppressed. Putting such kind of slabs into 
a series could create a macroscopic system with strongly enhanced pairing properties. More 
studies of this kind are under way. 
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